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SUMMARY 

A mathematical model is developed by which the following satellite 
orientation and control problems may be resolved: (1) determining at- 

titude for maximum area of solar cells in sunlight; (2) generating slew- 
ing commands for a change in attitude; (3) computing star tracker 
gimbal angles for maintaining proper orientation; and (4) determining 
when guide stars are occulted by the earth, sun, and moon. 
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MATHEMATICAL ANALYSIS FOR THE ORIENTATION AND 
CONTROL OF THE ORBITING ASTRONOMICAL 
OBSERVATORY SATELLITE 


by 

Paul B. Davenport 
Goddard Space Flight Center 


INTRODUCTION 

Sometime in 1964 the NASA expects to launch the first Orbiting Astronomical Observatory (OAO). 
The OAO, consisting of the spacecraft and the observatory (experimenter's package and equipment), 
is expected to have a lifetime of one year. The spacecraft, orbiting at an altitude of approximately 
500 miles, will carry highly refined equipment to conduct astronomical experiments free from the 
disturbing influence of the earth's atmosphere. 

The spacecraft will be capable of directing the experimenter's equipment anywhere in space and 
of maintaining this direction with a high degree of accuracy. This will be accomplished by a complex 
stabilization and control system whose major elements are: (1) sun sensors, rate gyros, and a jet 
system to initially orient and stabilize the spacecraft; (2) coarse inertia wheels to reorient the space- 
craft, and fine inertia wheels to maintain the desired orientation; (3) a set of six star trackers to de- 
termine the orientation of the OAO; and (4) a magnetic unloading system to keep the momentum in the 
fine inertia wheels from building up excessively. 

To place the OAO in the proper orientation in space, the precise slewing angles needed to accom- 
plish this feat must be known, as well as the gimbal angles of the star trackers when they are locked 
on their respective guide stars. In addition to these problems various physical constraints such as 
star occultation by various bodies including the earth, sun, and moon must be considered. There are 
also many constraints due to the spacecraft itself, which include such things as obtaining maximum 
power by proper positioning of the solar paddles, keeping the optical axis away from the sun’s dam- 
aging rays, and restricting the gimbal angles within their limits. This report describes the develop- 
ment of a mathematical model to aid in the solution of these problems. 


COORDINATE SYSTEMS 

The relative positions of the various bodies in space will be defined by giving their directions 
with respect to a fixed right-handed rectangular coordinate system with axes u, v, w and origin at the 
earth's center of gravity. The u-v plane will be the equatorial plane at some epoch T, the 
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positive u axis pointed toward the vernal equinox and the positive w axis toward the true north pole at 
time T. 

The position of a star will be denoted by its right ascension a and declination 8. In most applica- 
tions the direction cosines a, b, c of the star will be needed; these are given by 

a = cos 8 cos a , 

b = cos 8 sin a , 


c = sin 8 . 

Since the earth revolves about the sun, the position of the stars is not exactly fixed in this co- 
ordinate system; however, the maximum stellar parallax due to this motion is about 0.8 second of 
arc for the nearest star. 

The positions of the sun, moon, and OAO will be given at any time t (measured from the epoch T) 
by their respective orbital elements, which are: 

n = mean angular motion in the plane of the orbit 

e = the eccentricity of the orbit 

n = the right ascension of the ascending node; 0 5 ft < 2tt 
i = the inclination of orbital plane to the u-v plane; 0 1 i 5 tt 
a) = argument of perigee; 0 < a> < 2jt 
M 0 = the mean anomaly at time T 

Since no orbit is truly Keplerian, first and higher derivatives of the various elements may be 
necessary according to the accuracy required. It may also be necessary to update these elements 
from time to time. 

A secondary coordinate system with axes x, y, z is defined with respect to the OAO. The origin 
of the system is at the geometric center of the OAO with the optical axis coincident with the x axis. 
The positive slewing motions of yaw, pitch, and roll are defined to be clockwise rotations, as seen 
from the origin, about the positive z, y, and x axes respectively. 

The attitude of the OAO in space is described by specifying the right ascension a and declination 8 
of the x axis and the angle /3, which is the angle the y axis makes with the u-v plane measured in 
the y-z plane (positive direction toward the positive z axis). 

Because of the great distances involved, the angular coordinates of the stars and the sun will be 
considered to be the same whether the origin of the coordinate system is at the center of the earth or 
the spacecraft. In the case of the sun, placing the center at the OAO instead of the earth produces a 
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maximum error of less than 10 seconds of arc because of parallax (for a 500 mile circular orbit). 
For the stars this error is completely negligible. However, for the moon the error can become as 
great as 1 degree. 

For each star tracker on the OAO we define a coordinate system with axes x., y., z. such that 
the outer gimbal angle o is in the x. - y. plane (positive angle measured from positive x. towards 
positive y t ) and the inner gimbal angle ji is the angle from the x. - y. plane (positive towards nega- 
tive z.). Thus the direction cosines of a star (with gimbal angles a and in this coordinate system 
are given by 


a = cos /x cos cr , 


b = cos ju s in cr , 


c = -sin /x . 


TRANSFORMATIONS 


The mathematical approach here is based on the algebra of rotations (matrix algebra); but, since 
we are concerned only with directions, the rectangular coordinates of a vector will be identical to its 
direction cosines. Thus, from the preceding definitions, if the attitude of the OAO is given by a, S,/3 
and if (u 0> v 0 , w 0 ) is any vector in the u, v, w system, then the coordinates x 0 , y Q( z 0 of this vector 
in the OAO system are given by 


where 



In addition, if the spacecraft is yawed, pitched, or rolled by <9, or respectively, the coordinates 
x,, y x , zj of the vector in this new coordinate system are given by 
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where 



or 





'10 0 
0 cos <f> sin d> 

0 -sin<£ cos <pi 


In order that all angles y may be defined uniquely, we use the function 


y - tan 1 (a/b) (1) 

with the understanding that the sign of a is the same as that of sin y and the sign of b is the same as 
that of cos y. Thus the signs of a and b in Equation 1 determine the proper quadrant, and the inverse 
tangent of a divided by b determines the proper value. For all other inverse functions the principle 
value shall be taken. The range of all angles defined above is as follows: 


Right ascension of x axis a 

Declination of x axis S 

Roll of y axis fi 

Yaw of OAO 

Pitch of OAO 6 

Roll of OAO <t> 


0 < a < 277 

- 77/2 £ S £ 77/2 

-■77 £ fi < 77 

- 77 £ \p < 77 

-77 £ 6 < 77 

-77 £ <fi < 77 


DETERMINATION OF FINAL ROLL 

If a 2 and s 2 are the right ascension and declination of a new target star S 2 , then we wish to de- 
termine the roll p 2 which will orient the solar paddles such that they receive maximum sunlight when 
the optical axis points towards the star S 2 . This will be the case when the angle y between the sun- 
line and a normal line of the paddle plane is a minimum. Let u s , v s , w. be the direction cosines of 
the sunline in the u, v, w coordinate system. The direction cosines x s , y s , z s of the sunline in the 
x, y, z system (optical axis pointing toward S 2 ) are then found as follows: 
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The direction cosines of the sunline in the x, y, z system after a roll of p are given by 

X s * 

y s cos p + z s sin /3 , 
z cos /3 - y s sin /3 * 

If x y z are the direction cosines of a directed normal from one side of the paddle plane, the 
cosine of the angle y between this normal and the sunline is expressed by 

f (P) = cosy ( 2 ) 

= x p x s + y p (y s cos p + Z s sin P) + Zp (z s cos - y s sin /3) . 

The requirement that y be a minimum implies that cos y be a maximum; that is, 

df _ _dy_ 

d/5 " - sin y 


= y p (z s cos/3 - y s sin/3) - z p (z s sin/3 + y s cos^) 

= (y P z s _ z P y a ) cos P ~ (y P y* + z P z ») sin ^ 

= o , 


or 


There are two possible values of J3 that will satisfy Equation 3: fi 0 and v However, the value that 
maximizes Equation 2 is the desired value of p\ denote this value by /3 0 . 

The above analysis has considered only one side of the paddle plane, but the other side is handled 
in exactly the same way with the normal whose direction cosines are -x pf -y p , - z p . This requires 
that the negative of Equation 2 be a maximum, which again leads to Equation 3. The proper solution 
this time will be the other value of p to satisfy Equation 3, namely, /3j. Thus, to determine the best 
value of /3, we evaluate f(/3 0 ) and — f (/5 X ) , and pick the value of p that gives the maximum of the two. 

GENERATION OF SLEWING COMMANDS 

If the OAO has an initial attitude of a t , , p x and it is desired to reorient in order to obtain an 

attitude of a 2 , $ 2 , p 2 , the slewing commands needed to accomplish this reorientation must be 
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determined. Since the OAO may be rotated about any one of three axes, there are twelve possible 
slewing sequences. These slewing sequences are listed as follows: 


yaw - pitch - roll 
roll - pitch - yaw 
yaw - roll - pitch 
pitch - roll - yaw 
pitch - yaw - roll 
roll - yaw - pitch 


yaw - pitch - yaw 
roll - pitch - roll 
yaw - roll - yaw 
pitch - roll - pitch 
pitch - yaw - pitch 
roll - yaw 7 - roll 


The analysis for determining the amount of slewing required is similar regardless of the slew- 
ing sequence; therefore we shall refer to a general slewing sequence of i-j-k. The matrices of these 
i-j-k rotations will be denoted by T., T., and T k respectively. 

If v is any vector with coordinates given in the OAO system with attitude a 1$ Sj , /?,, the matrix to 
find the coordinates of v in the OAO system with attitude a 2 , § 2 , p 2 can be obtained as a product of 
six matrices: 


T “ 1 T “ 1 


where T" 1 indicates the inverse of T. Likewise, if a slewing sequence of i-j-k is given when the OAO 
has an attitude of a p , p t , which causes the spacecraft to have a final attitude of a 2 , S 2 , the 
matrix of the transformation from the initial attitude to the final is found from the matrix product 
T k T. T. . Therefore T., T. , and T k must satisfy the matrix equation: 


T.T . T . 

k J i 


= T* T x T 


“ 1 T “ 1 T “ 1 
*1 Ts l T *1 


(4) 


The right-hand side of Equation 4 is a 3 x 3 matrix that can be determined from a,, p ,, a 2 , s 2 , p 2 . 
Therefore this matrix is independent of the slewing sequence. This right-hand matrix shall be denoted 
as C with elements C. , Thus for each slewing sequence the left-hand side of Equation 4 can be com- 
pared with the C matrix to determine the amount of slewing. Several examples are given below: 

SEQUENCE: YAW- P ITCH-ROLL 


T * T e T 4, 


( cos 9 cos p , 

sin sin 6 cos 4> - cos 9 sin <p, 
cos 4> s in 9 cos p + sin^sin 


cos 9 sin 4 s * 

sin 4> sin 9 sin \p + cos$cos</>, 
cos <p> sin 9 sin pi - sin^cos^, 


- sin 9 
s in <f> cos 9 
cos d> cos 9 j 


and comparing with the c matrix implies that 


sin 9 - ~C 13 , sin^ 

os e = ±yt, 2 , + c , 2 2 = ± + c 3 2 3 , cos^ 


Cj 2 / cos & » sin? 
Cj x /\ cos 9 y cos 4> 


c 23/ /c°s e , 

c 33 /c°s e . 
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SEQUEN C E: ROLL- YAW- PITCH 

( cos 6 cos ip , cos 8 sini/i cos <f> + sin 8 sin $ , cos 8 sin \p sin cf> - s in 0 cos 4> 

- sin ./- , cos $ cos 4> , cos «/» sin 0 

sin 8 cos \p , sin 0 sin <// cos 0 - cos0sin</> , sin 0 sin i/> sin 4> + cos 8 cos 4> 

and comparing this with the C matrix implies that 

sin\p = -C 2 j , sin 0 = C 3 ,/cos </> , 

cos./- = ± -|/C l 2 1 + C/, = ± i/C 2 2 2 + C/ 3 , cos 8 = C, ,ycos tp , 

SEQUENCE: YAW- PITCH- YAW 

j cos ip 2 cos 8 cos 1 p i ”* sin^j sin \p 2 , cos \p 2 cos ^ sin i//j + sin cos 

T , T^T, = [ - ( sirn/f 2 cos ^cosi//j + sin^, cos^ 2 ) , - sin cos 0 sin 1/^ + cos co 

v 2 v v i \ 

y sin£cos<//, > sindsinxp^ 

and comparison to the C matrix gives 

sine = ±yc/, + C 3 2 2 = i/C, 2 , + c 2 2 3 , sin = C 32/ /sin 0, sinxf> 2 = 

cos 6 = C 33 , cos/., = C 3l/ 4ine, cos/- 2 = 

Thus for each slewing sequence we obtain two solutions corresponding to the plus and minus sign 
of the radical. Therefore there are actually twenty-four possible slewing commands. 

We could continue the above process for all twelve combinations of the matrix T k T. T. , compare 
each with the C matrix, and find that there is a definite pattern as to which elements of the c matrix 

to choose for a given slewing sequence. Let the numbers 1, 2, and 3 be used to represent roll, pitch, 

and yaw respectively (2-1-3 would indicate a pitch- roll- yaw sequence); then, for all sequences i-j-k 
where the same type of slew is not used more than once, we have the following relations: 


a l23 ~ CT 231 “ ^31 2 ~ “°132 

a 213 

= - CT 32I = 1 

sin j = 

a i jk ^ki • 


cos j = 

J, - 

+ c 2 

^kk ’ 

sin i = -a.. k C kJ /cos j , 

sink ^ 

- CT i ik C ii/ cos i 


cos i = C kk ^cos j , cos k _ C ii ^cos j 


i//j , - cos \p 2 sin 6 

s i p x , sin \p 2 sin 6 

. cos 8 


C 23 /sin 8 


' C j 3 /sin 8 . 


sin 4> - C 2i/ /co? i p , 
cos <p - C 22/ /cos i/> . 
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In the degenerate case where C k . = C k1t = C 4i = C.. = 0, the value of j is either plus or minus 
90 degrees. In this case i and k must satisfy the following: 


sin (k ± cr ijk i) = J , 

cos (k ± a. Jk i) = T^ iik C ik • 

The upper signs are taken when j is +90 degrees, and the lower sign when j is -90 degrees. 

Example 1 — Assume the matrix C is given. Determine the slewing angles for a pitch- roll-yaw 
sequence. 

By definition i is 2, j becomes 1, and k is 3. Here a is - 1; therefore, 

sin - -C 32 , 


cos 4> = ± yfc 2 2 2 + Cj 2 2 


sin 9 

'= C 3l/ OS< ^ ’ 

sin \p 

~ C 12 / 

/ cos <p 

cos 0 

= 3 S <£> , 

cos \p 

- r 
“ ^2 2 / 

/cos <p 


If the slewing sequence includes the same type of slew twice (of the form i t - j - i a ), and k is the 
slew not used, then the slewing angles are defined by 


sinj = ±|/c i 2 J + C.\ = ±|/C J. + C k 2 , 


j =c n , 


tij = C . j^/sin j , sin i 2 = C ^sin j , 

= a. .C., /sin j , cosi 0 = -cr..C../sin 

j i i k/ J 2 j i k i/ 


^12 “ ^23 “ 1 = 3 “ _ Cr 21 “ °"32 ~ 1 ’ 

When Cj . = c ik = C j; = c ki = o, the angle j is either 0 or 180 degrees, depending on whether C. 
is plus or minus. In this case the angles ij and i 2 must satisfy the relations 
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sin(i 2 i ij) = cr n C kj , 
c °s(i 2 ± i,) = ± C kk , 


where the plus sign is taken if c. . is positive and the minus sign if C i . is negative. 
Example 2 — Determine the slewing angles for a yaw-pitch-yaw sequence. 

In this case i becomes 3, j is 2, k equals 1, and a is +1: 


sin 6 

- + i 

/ C 3 2 2 + C /. 

- ± 

Jr 2 + C 2 

r^23 c i3 * 



cos 6 ~ 

: C 33 ■ 


sin0j - 

C 3 2/ 

/s in 6 , 

sin 0 2 

= C 2 3 / Sin0 

cos \p j = 

C 3 V 

/sin 0 , 

cos 4*2 

= ~ C 1 3 / Sin0 


DETERMINATION OF GIMBAL ANGLES 

The determination of the star tracker gimbal angles will depend on the physical mounting of each 
star tracker with respect to the OAO’s coordinate system. Let T. ( i = 1,- • •, 6) be the transformation 
that determines the coordinates of a vector in the OAO system from the coordinates of the vector in 
the i th star tracker system. Thus, if and /x. are the outer and inner gimbal angles respectively 
of the i th star tracker, the direction cosines in this star tracker system are given by cos /*. cos cr., 
cos /*. sin a. , -sin n . and the direction cosines in the OAO system are obtained from the following 
relation: 



Thus, if the OAO has an attitude defined by a 2 , p 2 and the i th star tracker is locked on a star with 
right ascension a. and declination S., the direction cosines x. , , 2 . of the star with respect to this 

star tracker system can be determined from 



/cos 8 . cos a,\ 


T. 1 T« Ts T cosS. sina, 
i 0 2 S 2 a \ 

\ 

s in S . 
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The gimbal angles a. and are then found from the following expressions: 


cr. = tan 1 (y , 

Mi = tan~' (-Ziy/^Xj 2 + y ; 2 ) - 


If the OAO is slewed with a j-k-f slewing sequence and gimbal angles o-. ^ /z. ? after the slewing 
are desired as a function of the gimbal angles cr. , before the slewing, the following formulas may 
be applied: 



( COS U. COS CT \ 

A 

cos^i, sinff i, 

- sin Mj j / 

ct i 2 -- tan_1 (y iy Ai) , 


Mi 2 = tan ’ 1 (-z,/V x 1 2 + y i 2 ) • 


SATELLITE CONSTRAINTS 

Of the twenty-four possible slewing sequences there may be several that are not applicable be- 
cause of restrictions imposed by the spacecraft itself. These restrictions may be investigated by the 
same techniques employed in the earlier sections. 

One such restriction is that the sun shade may not protect the experimenter's equipment from 
the sun's rays if the optical x-axis is within 45 degrees of the sun. The first OAO also will have an 
experiment at the opposite end of the optical axis; thus in this case the minus x-axis must also be 
kept 45 degrees from the sun. 

Although the sun shades are designed to shut if either experimental axis comes within the pro- 
hibited area of the sun, it is desirable to avoid a slewing sequence that would require such action. The 
avoidance of such slews will eliminate damage to the experimental equipment even if the sun shade 
fails to work properly. 

Another satellite constraint is that at least two star trackers must not exceed their gimbal limits 
during the entire slewing sequence. If this is not possible, new star assignments must be made at 
intermediate slews. 
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To determine whether a target star with right ascension a and declination « lies within the pro- 
hibited area of sun, we merely determine the angle from the star to the sun. The cosine of this angle 
y is given by 


cos y - x cos S cos a + y cos 8 sin a + z sin S , 

s ■' s s’ 

where x s , y g , z s are the direction cosines of the sun. Thus, if y is less than 45 degrees, this star may 
not be viewed. In the case of the double-ended OAO, y must lie between 45 and 135 degrees before 
viewing is permissible. 

To check whether a slewing sequence will cause the optical axis to come within 45 degrees of the 
sun, the angle y between the sun and optical axis can be written as a function of the slewing angle. 
From this functional relation determine the domain of the slewing angle that makes y = 45 degrees or 
less. If the desired slew lies within this domain, the slew is prohibited. If a , § , p define the attitude 
of the OAO before the slewing sequence begins and if u s , v s , w s are the direction cosines of the sun in 
the inertial system, the direction cosines of the sun in the OAO system x s , y s , z s are given by 



After a slew of X, where k may be either a yaw, pitch, or roll, the direction cosines of the sun 
x', y B '. z s after the slew are defined as 


X 1 
s 



/ 

z 

s 



The cosine of the angle between the sun and optical axis as a function of \is then 


x cos k + b s i n k , 

s 1 


where b - y s if k indicates a yaw and b = -z s if k is a pitch. A roll slew need not be considered, 
since a roll does not affect the angle between the sun and optical axis. 

If k 0 is the desired slew, the slew will be allowable if 


x s cos k + b sin \ < cos 45 ° 


(5) 
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for all K between 0 and \ Q . Equation 5 may be written as 


x 5 b 

— rns X + — Sin A 

ysT 2 + b 2 yV + b 2 


< 


cos 45 ° 

y^7 + b 2 


or 


cos 45 ° 

cos e cos A. + sine sin A = cos (e-A) < — 

yx s 2 + b 2 


where 


€ 


= tan' 1 ( b /* s ) 


- 'TT < e <77 


If 


cos 45 ° 

— > 1 

y/x s 2 + b 2 


the slew is always permissible. If 


cos 45 
/x s 2 + b 2 


the slew is permissible only if none of the following angles lie between 0 and A 0 


6 - 7 ] , 

€+ 7] , 

e - 77 + , 

e + 77 - , 


where 



The next slew in the sequence may be checked in exactly the same way after replacing x s , y s , z s 

b Yv y:> <• 

The gimbal angles after each slew of a sequence may be obtained in the same fashion that the 
final gimbal angles are determined. If cr. and /!. are the outer and inner gimbal angles of the i th 
star tracker and a slewing sequence of Z-m- n is to be performed, the gimbal angles after the j th 
slew are given by 


where 



OCCULTATION 


The knowledge of stellar occultation is necessary for several important reasons: First, the at- 
titude of the spacecraft cannot be maintained if less than two star trackers are tracking at any time; 
second, the occultation of the target star during an experiment would result in wasted time. 

The three bodies that cause stellar occultation are the earth, sun, and moon. Occultation due to 
the sun and moon will be less frequent than that due to the earth. Because of this, the time of oc- 
cultation due to the sun and moon can be determined in a similar manner. In both cases the origin of 
the coordinate system is assumed to be centered at the OAO. This introduces an error due to paral- 
lax of about 10 seconds of arc for the sun and about 1 degree for the moon. 
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Let i and ft be the inclination and right ascension respectively of either the moon’s or sun’s 
orbit, and let a. and S. be the right ascension and declination of the i th star. The direction cosines 
a, b, c of this star in a coordinate system whose x-y plane is the plane of the orbit are given by 




COS ft 

- s in ft 


sin ft 
cos ft 


0 




The cosine of the angle y between the star and the sun or moon is then 


cosy = a cos + + b sin (a j+/x) r 


where ^ is the argument of perigee and p is the true anomaly. The general requirement is to de- 
termine when the angle y will be less than some fixed angle \ (45 degrees for the sun, one-half the 
angle subtended by the moon plus errors for the moon). When y just equals the corresponding time 
t is the time of immersion or emersion for that star. This requirement of equality may be expressed 
as 


a cos {<xi + fj.) + b s in [co + jjl) ~ cos K 


(6) 


or 


cos e cos ioo + fj.) + sine sin {co^/jl) 


cos \ 

+ b 2 


(7) 


where 


e = tan -1 (b/a) , 0<e<277 

Equation 7 may also be written as 


cos ( e - co ~ ft) 


cos \ 

+ b 2 


so that 


M = € ** CO* T) 


where 


V 


cos \ 

+ b 2 


+ 2krr 


k - 0,1,2,'’ - 
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Thus each value of v determines the value of m at an immersion or emersion of the i th star. This 
value of fj. will correspond to an immersion if 


sin(e-cy“/x)>0 , 


and to an emersion if 


sin { e - co- (i) <0 


Once fi is known, the corresponding value of time can be obtained by the following relations: 


E 


tan 


fl 1 - e 2 sin jj. \ 

, 9 T-^-l 

\ e + cos fi J 

+ 277 L 277 J 


M = E-esinE 


t 



( 8 ) 


where e is the eccentricity of the orbit, n the mean motion of the body in the plane of the orbit, M 0 the 
mean anomaly at epoch, and t the time from epoch. 

Although occultation due to the earth is more frequent and troublesome than that due to the sun 
and moon, it can be handled in a similar manner. The angle \ (immersion or emersion occurs when 
the angle between OAO and the star equals k) in this case is not a constant but a function of the OAO's 
range. Thus immersion and emersion can be determined analytically only when the spacecraft is in 
a circular orbit. For small eccentricity, however, the circular solution using the mean range should 
give sufficient accuracy. 

Because of the great distances of the stars we may assume that the line from the origin to a star 
is parallel to the line from the OAO to the star. Then, with the additional assumption of a spherical 
earth, the cosine of k (see sketch) may be determined by simple trigonometry. 



To star 


[x] indicates the greatest integer less than or equal to x. 



Thus 


cos \ 


-yp^r 


r 


where the unit of distance is the earth’s radius. Hence, if r is constant, occultation by the earth may 
be handled in the same manner as that of the sun and moon, that is, by Equations 6 through 8, where 
the elements are those of the OAO’s and where the x-y plane of the coordinate system is the OAO’s 
orbital plane. The test for immersion or emersion in this case is the reverse of that given for sun- 
moon occultation. 

CONCLUDING REMARKS 

The analysis and mathematical models contained in this report are intended to be quite general. 
In many cases the formulas may be simplified if accuracy requirements warrant it. In other cases 
a different coordinate system will simplify some models; for example, the ecliptic system would re- 
duce any model involving the sun’s coordinates. Thus the formulas contained herein are not de- 
pendent on any particular coordinate system, and simplifications may be made by simply omitting 
terms in various expressions. 
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